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^! Abstract 

We calculate neutrino emissivities from the decay and scattering of Goldstone bosons 
in the color-flavor-locked (CFL) phase of quarks at high baryon density. Interactions in 

> 

the CFL phase are described by an effective low-energy theory. For temperatures in the 

00 ' 

tens of keV range, relevant to the long-term cooling of neutron stars, the emissivities 

m 

■ involving Goldstone bosons dominate over those involving quarks, because gaps in the 
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CFL phase are ~ 100 MeV while the masses of Goldstone modes are on the order of 
10 MeV. For the same reason, the specific heat of the CFL phase is also dominated by 
the Goldstone modes. Notwithstanding this, both the emissivity and the specific heat 
from the massive modes remain rather small, because of their extremely small number 



^ ■ densities. The values of the emissivity and the specific heat imply that the timescale 

for the cooling of the CFL core is ~ 10 26 y, which makes the CFL phase invisible as 
the exterior layers of normal matter surrounding the core will continue to cool through 
significantly more rapid processes. If the CFL phase appears during the evolution of 
a proto-neutron star, neutrino interactions with Goldstone bosons are expected to be 
significantly more important since temperatures are high enough (~ 20 — 40 MeV) to 
admit large number densities of Goldstone modes. 

PACS: 26.60,+c, 95.30.Cq, 97.60.Jd 
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I. INTRODUCTION 



Neutrino emission and interactions in matter at high baryon density play crucial roles in astro- 
physical phenomena such as core collapse supernovae and neutron stars [|lj] . Neutrinos drive supernova 
dynamics from beginning to end: they become trapped within the star's core early in the collapse, 
forming a vast energetic reservoir, and their eventual emission from the proto-neutron star is prodigious 
enough — containing nearly all the energy (~ 10 53 ergs) released in the explosion — to dramatically con- 
trol subsequent events [Q]. The //-luminosity and the time scale over which us remain observable from 
a proto-neutron star are also governed by charged and neutral current interactions involving matter 
at high baryon density ||. The long-term cooling of a neutron star, up to a million years of age, is 
controlled by i/-emissivity and the specific heat of the densest parts of the star; during this time the star 
is observable through photon emissions, which may allow us to determine the star's mass, radius, and 
internal constitution || . The tabulation of surface temperatures and ages of neutron stars is currently 
one of the primary goals of X-ray astronomy. 

Collins and Perry || noted that the superdense matter in neutron star cores might consist of weakly 
interacting quarks rather than of hadrons, due to the asymptotic freedom of QCD. Asymptotic freedom 
implies that at very high baryon density, for which the baryon chemical potential /i_b 3> Aqcd, QCD is 
amenable to perturbative techniques. However, recent studies have shown that the naive ground state 
of the system, a Fermi liquid of weakly interacting quarks, is unstable with respect to the formation 
of diquark condensates |6|-|9| . Attractive interactions induce the quarks to pair which results in a gap 
in their excitation spectrum that may be as large as A ~ 100 MeV [T^JTl|. The quark Cooper pairs 
are anti-symmetric in color. If only two flavors participate, the pair wave functions have the structure 
e abc u b d c , where a, b, and c are color indices. This implies that color SU (3) is broken to SU (2) by a 
Higgs mechanism and that one type of quark and three types of gluons remain massless. These degrees 
of freedom will dominate the low energy excitations of the systemF]. If the baryon chemical potential 



If the temperature is very small, T < 1 keV, the ungapped quark may acquire a gap and the gluons will be 
confined into glueballs. 
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Hb is very large, we expect that all three light flavors participate. In this case, the wave function of the 
Cooper pairs is expected to be color-flavor locked fl2|-|T5[|. Color-flavor locking implies that all quarks 



and gluons acquire a gap. 

Because the critical temperature associated with color superconductivity is large, we expect that 
quark matter existing in the core of a neutron star will be in a superconducting state. This leads to 
the question whether one can infer the presence of a color superconducting core from astronomical 
observation, and whether such observations can distinguish among the different color superconducting 
phases. 

Color superconductivity is a Fermi surface phenomenon and its effect on the equation of state is 
expected to be small. The impact of a diquark condensed phase on the long-term cooling of neutron 
stars was studied in |T6|Ji7|] . These works explored the consequences of neutrino emission processes 
that involved gapped quarks. The cooling of a pure quark star, considered in [16[], implied that cooling 



would occur very fast, indeed too fast to be consistent with the existing X-ray data. In [jD]], in which a 



hybrid star with a mixed phase of hadrons and quarks was considered, large quark gaps rendered quark 
matter invisible and vanishing quark gaps led to cooling behaviors which were indistinguishable from 
those of normal stars. A possible connection between the occurrence of quark matter in neutron stars 
and neutron star glitches was proposed in ]T8[] . Neutrino diffusion in the two-flavor superconducting 
(2SC) condensed phase was investigated in 



In this paper, we focus our attention on the long-term cooling stage, which begins when the star's 
temperature has dropped to a few tenths of an MeV (~ 10 9 K or T 9 = 1) and when most of its 
lepton content has been lost due to neutrino radiation during the so-called proto-neutron star (PNS) 
evolution that lasts about a minute or so. The extent of a quark phase is maximized when the trapped 
neutrinos have left the star p0| . Thereafter, the star continues to lose its energy by radiating low- 
energy neutrinos until it enters a photon cooling epoch and its surface temperature may be estimated 
by the detection of X-rays. Coupled with an independent estimate of the star's age, the cooling history 
of the star may be inferred. 

In a phase in which all quarks are gapped, the neutrino emissivity from the direct Urea processes 
fi + £ — > /2 + ui, fi — > fi + £ + V&-, where f\ and /2 are quarks and t is either an electron or 



muon, are strongly suppressed due to the energy cost (~ A) involved in breaking the pair. (The 
modified Urea process, which contains an additional quark in both the entrance and exit channels, is 
additionally suppressed.) However, there are other mechanisms for neutrino emission from this phase. 



One such process is the pair breaking and recombination of quarks ||21|| . This process can dominate 
over conventional cooling modes just below the transition temperature T c ~ A/2. We demonstrate 
here that as the temperature falls below T c , neutrino emission will be dominated by decay of massless 
or almost massless collective modes, hereafter denoted by <j), that are associated with the breakdown 
of a global symmetry. The low energy excitations of CFL quark matter are an exactly massless 
mode associated with the breakdown of the U(l) of baryon number and a very light octet and singlet 
of Nambu-Goldstone bosons associated with the breaking of the exact SU(3)l x SU(3)r chiral and 
approximate U(1)a axial symmetry. 

The purpose of this work is to provide benchmark estimates of the neutrino emissivity from the 
electroweak decays of Goldstone modes. In Sec. II, we employ an effective theory of CFL Goldstone 
modes coupled to electroweak gauge fields. Specifically, we calculate the neutrino emissivity from the 
processes 

1. TT ± l ± + V l , K± -> l ± + U l , 

2. 7r°(?7, 7/) — ► v + v , and 

3. + 0^0 + Z/ + Z/ , 

and compare our results to those from previously studied v— emission processes. In Sec. Ill, we 
calculate the photon emissivity from the pion electromagnetic radiative two-photon decay, tt° ^7 + 7. 
This is followed by a discussion of photon- Goldstone boson interactions in the CFL phase. The cooling 
behavior of a neutron star is also governed by the specific heat of the ambient matter. We therefore 
assess the contribution of the Goldstone modes in Sec. IV. Our principal findings are summarized in 
Sec. V. Utilizing our results for the emissivity and specific heat, we estimate the time scale over which 
a CFL phase core cools in the concluding Sec. VI. Here, we also comment on the expected role of 
neutrino interactions with Goldstone bosons in the evolution of a PNS, in which temperatures in the 
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range of 20-40 MeV are encountered. 



II. ELECTROWEAK INTERACTIONS OF GOLDSTONE MODES IN THE CFL PHASE 

The low energy excitations of the CFL phase are two singlet modes associated with U(1)b and U(1)a 
symmetry and an octet of Goldstone modes associated with chiral symmetry breaking. The octet is 
described by a low energy theory that bears strong resemblance to chiral perturbation theory p2Vp4 



The Goldstone modes are parametrized by a 3x3 unitary matrix U which is a color singlet, transforming 
under SU(3)l x SU(3)r as U — > gLUg R , where U is related to axial-like fluctuations of the left and 
right handed diquark fields: 

L ai ~ e abc e l3k {q b iqt)\ R ai ~ e a6 V ifc (g^ fe }*, and U — LR) . (1) 

The fields L and R carry color (ijk) and flavor (abc), and transform under gf C SU(3)f and g c C 
SU(3)c (where / denotes left or right handed flavor) as 

L -> giLgl and R -> g R Rgl , (2) 

respectively. The low energy effective theory is governed by the Lagrangian 

f 2 

■J IT 



d Ud U^ - vidiUdirf 



+ -.., (3) 



where v n is the Goldstone boson velocity and we have suppressed mass terms and higher derivative 
terms. The chiral field U can be related to the octet meson field by U — exp(i7r a A a // T ), where the 
SU(3) generators A a are normalized as tr[A a A fe ] = 25 ab . 

In the weak coupling limit, we have v\ — 1/3. As a consequence of the breaking of Lorentz 
invariance in matter, there are two pion decay constants, fx and f$, which characterize the coupling 
of the pion to the temporal and the spatial components of the axial- vector current. The effective 
Lagrangian in Eq. (|3|) implies that fx = f n an d fs — v nfn- This result is consistent with axial- vector 
current conservation in the limit m n — * 0. We observe that /W 2 — fsQ 2 — for an on-shell pion. 

The extension to include electroweak interactions is achieved by minimal coupling ||26||: 
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D^L = d^L - ieA^QL - -j= 



L 



<g 



2cos# 



w 



a 3 - 2Qsin 2 6 w 



(4) 



D^R = d^R - ieA^QR + 



—^—Z° 
2cos0 w M 



2Qsm 2 9 w 



R. 



Here, A is the electromagnetic fieldQ, which couples with strength Q = diag(2/3, —1/3, —1/3) to the 
u, d, and s quarks, W and Z are the SU(2) electroweak gauge fields, g is the weak coupling constant, 
and 9w is the Weinberg angle. The a matrices are the Pauli matrices of SU (2) embedded into flavor 
SU(3): 



1 \ 1 0^ 



a 








^0 0^ 



and a 3 
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1 
0-10 
0-1 



(5) 



^ = 1 


In the charged current coupling, we have only taken into account the first two flavors |26fl . We 
will return later to comment on charged current decays of kaons. We observe that the CFL pions, 
kaons, and etas have charged and neutral current interactions that are very similar to those of ordinary 
mesons. An important difference, however, is the fact that CFL mesons are very light, itlgb — 10 



MeV |24| ; see p7| , |28[| for a recent discussion. This implies that the pion and kaon are lighter than the 
muon and dominantly decay into electrons and neutrinos. Another difference, as compared to the zero 
density case, is the fact that the Goldstone boson dispersion relation is modified from its vacuum form: 



e p = Mf + y m GB + Kp 2 ■ 



(6) 



Here, fi e ff is an effective chemical potential that depends on the quark masses |29|| . For pions, fi e jf ~ 0, 
but for kaons /i e // is sizeable and may have consequences. The modified dispersion relation Eq. (0) 
leads to a number of unusual effects. One of them, the fact that a very fast pion becomes absolutely 
stable, was already discussed in |J{J. Here, we will encounter a second one, the presence of the "helicity 
forbidden" decay 7r° — > uu. 



2 The A field may be re-expressed in terms of the physical (massive) gluon and physical (massless) photon 
that arise as a consequence of photon-gluon mixing in the CFL phase. 
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A. Neutrino Emissivity from n^ip) — ► I (pi) + vi{P2) 



Substituting Eq. (|J) into the effective Lagrangian, we obtain the relevant interaction term for 
charged pion decay: 

Cw, = 9 ~y (w ~d\ + + W Q + d\-^J - (f T - fs, - •) . (7) 
The pion decay constant has been estimated using various methods in the literature fl2"4|j25| . In 



Ref. [24], it was computed by matching a chiral effective theory with a microscopic theory of quasipar- 



ticles and holes near the Fermi surface. Explicitly, 

2 _ 21-81n2 p? 
JlT 18 2tt 2 ' 



In the following, we will outline the calculation of the v— emissivity for the 7r decay. The 7r + decay 
gives exactly the same result. The T-matrix element for the weak decay process is 

(r(pM(p 2 )\T\7r-(p)) = (2n) 4 6\p - Pl - p 2 )G F u l (p 1 )( 1 °p f T - 7 • p/s)(l - 7 5 )^(p 2 ) , (9) 

where we have made the identification Gp/Vz = g 2 /(8My V ). Using the equation of motion for the 
leptonic fields, the neutrino emissivity from this process is given by 

£ - = / wk: n ^ I (4fc Mpi) / i2$k; Mm >> ~ pi - P2) £ |M|2 - (10) 

Y,\M\ 2 = 8G 2 F MelV • p 2 + \x\ 2 KK^ P2 + Pi • Pa) + mlipx ■ p 2 - 2u p u P2 )} 

spin 

+ 2|^||x|m e (a; p a;p 2 -pi-p*)), (11) 



where u p = \jrn l GB + v^p 2 , £, = —im e f T and x = ~ *(/t — fs)- Here, m e is the electron mass, ub{uj p ) 
denotes the Bose occupation factor for the pion, and np{pi) = 1 — np(pi) is the Pauli blocking factor 
for the outgoing electron. For the temperatures of relevance to long-term cooling, Pauli blocking of 
the outgoing neutrino may be neglected. For bulk CFL matter at T = 0, the electron chemical /i e — > 



p9| , |3Tl| . At non-zero temperature, there is a small electron chemical potential due to the fact that the 
energies of positively and negatively charged Goldstone modes are not the same. For T < 1 MeV, we 
find /i e <T and the effect on neutrino emission is negligible. In proto-neutron stars, CFL phases with 



charged meson condensates can also exist p2| . Since our focus here is on long-term cooling, we do 
not consider these possibilities. Finally, we expect that CFL matter inside a neutron star is in contact 



with a hadronic or quark phase that has a large electron chemical potential. As discussed in ||33|| , this 
will lead to a thin charged surface layer which shields the CFL phase. In the following, we will assume 
that /i e = in regions well separated from this layer. 



The squared matrix element in Eq. (|Tl|) consists of three terms. We will evaluate the emissivities 
from each in turn. The first term involves |£| 2 and has the Lorentz invariant structure familiar from 
its vacuum counterpart. This is the only term that survives in the limit fr = fs- We also observe that 
the emissivity is proportional to the electron mass squared, as expected. The momentum dependence 
is of the form p! ■ p 2 = u pl u P2 (l — viv 2 cos9) with v% = \p 1 \/u Pl , v 2 = \p2\/u P2 , and 9 = Z(pi,p 2 ). The 
angular integral gives 

J dQ p dQ Pl dQ P2 {l - ^ 2 cos#)(27r) 3 5 3 (p - Pl - p 2 ) (12) 
= 2w pl ^ P2 - (p 2 -p\ -p\) ©(p + pi -P2)®(p + P2 -Pi)6(pi +P2 ~p) , 

PPlP 2 UJ pi U) p2 \ J 

where p is now to be understood as |p| etc. The triangle inequalities imposed by momentum conser- 



vation restrict the range of the p 2 and p (or equivalently u p ) integrals as follows: 

uj p -p m 2 e u p + p m 2 e 



<P2< 



2(u p + p )-"- 2 2(u p - p) ' 



ml "-^ • (13) 



\ 1 

When these conditions are satisfied, the p\ integration can be done trivially using the energy-delta 
function. The expression for the emissivity then becomes 

A r*** j / ^( m l 2 K-i) 2\ r p ™ x j 1 

H = RA 3 2 / du p n B (u p )[ — -m e + — u ) . dp 2 p 2 — — , 14 

bA^v^Jmn \v* ' /Jp™ in 1 + exp [p 2 — ujpj/l 

where A = 8Gpf^ml and the limits co>™ aa: , p 2 nm , and p™ ax are defined by Eq. (|l"3"D. In obtain- 
ing Eq. (Q), we have utilized the pion dispersion relation u p = v 2 p 2 + m\ to change the variable of 
integration as pdp/u p = doj p /v 2 . We have also used the identity h F (x) = n F (—x). For long-term 



cooling, temperatures range from tens of keV to eV, much less than either the pion mass or elec- 
tron mass (~ 10 MeV and 0.51 MeV, respectively). In the following, we will therefore calculate the 
emissivity in the low temperature limit when m n /T and m e /T ^> 1. Unless the pion mass becomes 
very small, m n ~ m e , we can also neglect the electron mass in the integrand of Eq. (|T4|). Within the 
range of the u p integration, we always have p 2 < uj p . In the low temperature limit, m^/T ^> 1, we 
can replace the Pauli-blo eking factor x) by unity. In this limit, we can also replace the Bose 
distribution function for the pion by the Boltzmann distribution. Introducing the scaled variables 
x = ujp/T, y = P2/T, and if) = m n /T, Eq. ([14]) now reads 



AT 5 r^ma, (l-vl 



QAn 3 V 2 Jif, \v% \ v 



^ = ^2/. dxe ~ x \h- K / ( 15 ) 



Umax 



!/.■!' 1 !: 



where ipmax = 4>/ \f 1 — v%. The y integration yields x\Jx 2 — ip 2 /2v 7T with the result that 

AT 5 ftma* _ x XVX 2 - if) V , 2 2 , 2 \ 

^ = 7^K , dxe -^-U 2 -(l-vl)x 2 ) . (16) 



647r 3 ^ 2 
As ip is parametrically large, most of the contribution to the integral comes near the lower limit. The 
upper limit can be extended to infinity with little impact on the numerical value of the integral. Using 
the substitution x = ip cosh 9 (the limits on 9 now run from to 00), the definition of the modified 
Bessel function of the second kind K„ and their recursion relation, we finally obtain 



A 2 



1 /m,T\ 3 / 2 



1 2 ' (17) 



64tt 



where we have used the asymptotic form K u {z) = \Jn/2z e~ z valid for any v when z 3> 1. The factor 
A is proportional to the interaction strength and the factor in the square bracket is the number density 
n n of thermal non-degenerate pions. Note that the dispersion relation in Eq. (^J) introduces an extra 
factor of in the number density of pions relative to the case with v n = 1. 

The other two contributions to the emissivity in Eq. ( |TUD are evaluated along similar lines. The 
result for the mixed term (|£||x|) and the quadratic term (|x| 2 ) are 

and (18) 

respectively, with B = WGpfrifT — fs) m l and C = 16G^(/t — fs) 2 fn n T. We observe that the 
coefficient C is independent of the electron mass. The term involving C arises purely as a consequence 



of Lorentz symmetry breaking in dense matter. We also note that all emissivities appear to diverge as 
v n — ► 0, which is the case near a phase transition. This is an artifact of the approximation ip max — ► oo 
in Eq. (|T6|). At fixed T, all rates are finite as tv — ► 0. Combining Eqs. (|T7| ) and (I8"|), and using the 



perturbative result v 2 = 1/3, we obtain the neutrino emissivity from the electroweak decay of CFL 
pions as 

where Sf n = fr—fs- The rate from tt + decay is exactly the same. This result is valid when temperatures 
fall below an MeV. Numerically, the emissivity may be expressed as 

T 



e?r = 1.18 x 10 28 /4o m 2 10 (m 10 T 9 ) 3 / 2 e 



116m^Q 



1 + 1.14 



erg cm s , (20) 



where /iioo is the quark chemical potential in units of 100 MeV, mio is the pion mass in units of 10 
MeV, and T 9 = T/10 9 K. 

It is interesting to compare the neutrino emissivity from pion decay with the emissivity from the 
quark direct Urea process in a gapped superfluid, as well as the quark pair breaking and formation 



(PBF) process. From |34] , |35| , the emissivity from the quark direct Urea process is given by 



e qP ~ 8.8 x 10 26 a s f^Vv 3 ^ 6 e~ A ^ T erg cm" 3 s" 1 , (21) 
V no / 

where a s = g 2 /An is the strong coupling constant, tib is the baryon density, n = 0.16 fm -3 is the 
nuclear equilibrium density, and Y e = ub/uq is the electron concentration. For n# = 5no, a s ~ 1, 
Y e ~ 10" 6 , A ~ 100 MeV and m n ~ 10 MeV, we find that pion decay dominates the Urea emissivity 
for T < 5 x 10 10 K. The emissivity from the quark PBF, which is effective in rapidly cooling the star 
during its early hundreds of years, is given by |2l| 

1.4 x 10 20 N V T 7 9 F a q (— V 3 erg cm" 3 s" 1 , (22) 

V n / 



q ~ 



where N u is the number of neutrino flavors, F is a temperature dependent function of order one that 
vanishes exponentially with the gap, and a q is a flavor dependent numerical factor of order 0.1. In 
the vicinity of T c , there is no exponential suppression, but as T < T c , the function F approaches 
exp(— 2A/ ksT) so that the emissivity from PBF becomes very small. 

10 



B. Neutrino Emissivity from iT G [r],rf) — > v + v 

In the previous section, we noticed that the emissivity for 7r~ — ► e~ + v process in matter contains 
a term that is not proportional to vn? e . This implies that the neutral current decays 7r° — > ^ + z/ and 
77(77') — > 1/ + z/, which are forbidden by helicity selection rules in vacuum, can occur in matter. The 
reason these decays are forbidden in vacuum is that the two neutrinos have opposite chirality. Thus, 
the two neutrinos have total angular momentum one in the rest frame of the decaying meson, Hence, 
a scalar meson cannot decay into vv. In a dense medium, it is still true that a scalar meson which is 
at rest with respect to the medium cannot decay into vv pairs. However, boost invariance is lost. This 
implies that the wave function of a scalar meson which is moving with respect to the medium can have 
higher angular momentum admixtures when viewed in its rest frame. 

In the following, we calculate the emissivity due to the process tt° —>■ v + v. The interaction term 
relevant for this decay channel is given by 

*W = ( z o^°) ~ [h - fs, - i). (23) 



2 cos 6^ 

The corresponding interaction for the 77 decay has an extra factor of l/y/3 (l/v^6 for 77' decay) from 
the normalization of the 577(3) generators. For simplicity, we ignore tt° — 77 — 77' mixing |2^j28[| . We will 
outline the emissivity calculation for the ir° decay. The corresponding emissivity for the 77 is obtained 
by changing the overall factor and the mass. The neutrino pair emissivity from 7r° decay is given by 

£ - = / wkc aM I pfe / w?k}^ + ^' s ^ - pi - ^ £J M i 2 < 24 > 

£ \M\ 2 = W^xf^p^ + P, ■ P2) , 

spin 

where N v counts the number of neutrino flavors. The emissivity is proportional to fr — fs, so it is 
directly proportional to Lorentz symmetry breaking in matter. The integrals are performed in the limit 
m e — > and m^/T ^> 1 as before, and the emissivity from this process is found to be 

Nv —{G 2 F {8U) 2 ml)Tm,n„. (25) 



47T V\ 

Numerically, the emissivity (for one neutrino flavor) is 

e vP = 4.52 x 10 2 V? 00 m lo (™ioT 9 ) 5/2 e ~ erg cm" 3 s" 1 . (26) 



llCm 



11 



C. Emissivities from Kaon and Massless Goldstone Boson Decays 



In this section, we would like to provide a brief survey of other processes that contribute to the 
emissivity of CFL matter. Neutrino pair emission from plasmon decays in a stellar plasma, first 
proposed in |36| and subsequently investigated in detail in [^,^], can be the dominant energy loss 



mechanism for very hot and dense cores of red giant stars or white dwarfs. Blaschke et al. [fil| studied 
the analogous plasmon decay of the massive gluon in the CFL phase. They found that the rate is 
suppressed by exp(— m g /T), where m g ~ gfi is the effective gluon mass. This implies that the plasmon 
rate is even more strongly suppressed than the quark Urea rate. 

The physical photon A in the CFL phase is a linear combination of the ordinary photon (A) and 
gluon (G) fields 

Afj, = A^ cos 9 CFL ~ Gj sin 9 C fl, tan 9 CFL = 2e/V3g , (27) 

where e and g are the electromagnetic and strong coupling constants. The massive gluon G 8 is the 
orthogonal combination. In an electron plasma, the photon is dressed by particle-hole excitations and 
acquires an effective mass. As a consequence, it can decay into vv pairs. In CFL matter at T = 0, the 
photon is dressed by both particle-hole and particle-particle excitations, and the dielectric constant is 
larger than one, but it does not acquire a mass |59|,[|(J. As a result, it cannot decay into vv pairs. At 
T 7^ 0, the photons acquire a small mass, mainly because of thermal e + e~ pairs, but the corresponding 
7 — > v + v rate is very small. 



In CFL matter, the usual spectrum of meson masses is partially inverted pi| , p9f . In particular, we 
expect that Goldstone bosons with positive strangeness, the K° and the K + , are lighter than non- 
strange Goldstone bosons such as pions and etas. The K + contributes to neutrino emission via the 
decay K + —>■ e + + v. The emissivity for this process can be computed along the same lines as the 
emissivity from charged pion decay. However, there are two differences. First, kaon decay is suppressed 
by the Cabbibo angle sin 2 (9c), where 9c ~ 15°. The second difference is related to the fact that the 
kaon dispersion relation is modified by the effective chemical potential term, see Eq. (|6]). We find 

ejt± = ^El ± n K , (l + 2 If ± fii) + ± f^V ) , (28) 

64vr ^ \f v E K ± J v^mj \ f w E K ± J 
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where Ak = 8 sin 2 {6c)G 2 F f^m 2 is the effective coupling, and E K ± = rriK =F faff with /i e // = m 2 /{2jj, q ) 
is the energy of a kaon at rest. Note that in the limit m e — > 0, 5f w — > 0, and — > 1, the emissivity 
remains finite because the terms containing /j, e ff breaks Lorentz invariance. In this case, 



sm\e c )G 2 F f 2 

K± = W" m K±n K ±Tn eff (29) 



For T ~ 10 9 K, this contribution to the emissivity would dominate over that of the first two terms 
in Eq. ( |28| ) that are proportional to ml . We note that the K + emissivity is not equal to the K~ 
emissivity This implies that, strictly speaking, it is not legitimate to ignore the effects of a non-zero 
electron chemical potential. We also note that, because of the smaller energy of a kaon compared to a 
pion, kaons will likely dominate the emissivity despite the Cabbibo suppression. Numerically, we find 

116E 10 

e K = 3.38 x 10 26 /4o E 2 W (m 10 T 9 ) 3 / 2 e~ — 

2, 



X 



l i 2 r ± Ms 1 + 2niK±T ( s f* ± ^ff 



erg cm s , (30) 



f n E K ± J v 2 K m\ \ f n E K ± j 

where m 10 and E w are the kaon mass and energy in units of 10 MeV. In Sec. [I FBI, we saw that in 
CFL matter the decay tt° — > v + v is allowed if the 7r° has non-zero momentum. The analogous decay 
K° — * v + v, however, remains strongly suppressed because it requires a second order weak s — > d 
transition. 

Since the CFL phase is characterized by unusually light Goldstone bosons, there is a possibility 
for Bose condensation to occur. Indeed, it has been argued that CFL matter at densities that can be 
achieved in neutron stars is likely to support a K° condensate p9|j41|j42]| . Furthermore, because of the 
presence of trapped neutrinos, CFL matter in a proto-neutron star may have charged pion or kaon 
condensates |[32|| . In neutron matter, pion condensation (if it occurs) leads to a substantial increase of 
the neutrino emissivity through processes like n + 7T + ^n + e + + v [f43| , |44f . In CFL matter, there are 
no ungapped fermions and the analogous process is suppressed by exp(— A/(fcsT)). 

If charged pions or kaons become lighter, then the emissivity from the decay process tt^ — > e ± + v 
will initially be enhanced. There is, however, no neutrino emission from the decay of a massless 
charged pion or kaon in a Bose condensed phase. The reason is that the decay of a massless boson 
with dispersion relation to = ||&| into pairs of leptons is kinematically forbidden. Indeed, we observe 
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that the emissivity for the process 7r =l= — > e ± + v is of the form m^ 2 exp(— m w / (fc^T)) and vanishes as 
— > 0. 

Instead, we have to consider neutrino Bremsstrahlung in Goldstone boson scattering, or Goldstone 
boson scattering followed by the decay of a Goldstone boson into vv or e ± z>. In the case of massive 
Goldstone bosons, these processes are suppressed by an additional power of exp(— m GB / {k B T)) with 
respect to the direct Goldstone boson decay process. For massless Goldstone bosons, however, neutrino 
emission is dominated by Goldstone boson scattering. 

In CFL matter, there is always at least one exactly massless Goldstone boson which is associated 
with the breaking of the U(1)b of baryon number. The U(1)b Goldstone boson contributes to the 
neutrino emissivity through the process + — > + z/ + z/. The coupling of to the Z-boson is 
given by (g/(12 cos(9w))Z fl fd lJ '(j) where / is the U(1)b decay constant p4|. The tctt and KK scattering 



amplitudes are fixed by the leading order 0(p 2 ) chiral Lagrangian. The 00 scattering amplitude, on 
the other hand, vanishes at leading order in the low energy expansion. The 00 scattering amplitude is 
on the order of A ~ p 4 /{f 2 fJ' 2 ) ||45|| . The emissivity will then scale as 

<°°~%r- pi) 

Including numerical factors associated with the scattering amplitude and the phase space integral, 

e GB ~ KT 11 x T 9 15 /ir 6 erg cm' 3 s" 1 . (32) 

This estimate is adequate for T <C A. For the range of temperatures considered in this work, this 
process is not likely to play an important role. Since there is no exponential suppression, it will, 
however, dominate neutrino emission at very late times. 

III. PHOTON EMISSIVITY AND PROPAGATION IN THE CFL PHASE 

So far, we have identified sources for neutrino emission from CFL quark matter. Another interesting 
issue is the photon emissivity from this phase. A first step towards computing the emissivity of 
photons from the CFL phase was taken in ||40|| , where the contributions from pair-correlated quarks 
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and qq annihilations were found to be comparable in magnitude to the rates from a hot hadron gas, at 
temperatures of several tens of Me Vs. At keV temperatures, however, the rates are vanishingly small 
due to the fact that quark gaps are of order 100 MeV. 

Another source of photons in CFL matter is the anomalous 7r° — > 7 + 7 decay, which proceeds via 
the non-abelian axial anomaly. The photon emissivity due to this process is given by 

1 rd?qn(u q ) f d 3 p{l + n{u p )) f d 3 k (1 + n(u k )) 44 
^|A/( Jr "^27)r = 24(p.t) 2 . 

spin 

The factor 1/2 accounts for the two identical bosons in the final state. The distribution of photons in 
the final state is Bose enhanced. The constant Ay is related to the coefficient of the electromagnetic 
anomaly for the isospin axial vector current. Anomalous electromagnetic processes in the CFL phase 
can be studied in analogy with the vacuum case |26|,|46j through the replacement e — > e = ecos#. We 
have 

« 3 = -^^F Q/3 F^, (34) 

which leads to the identification Ay = e 2 /47r 2 / 7r []4"7 |. Unlike in vacuum, the matrix element (p.k) 2 now 
depends on the pion momentum: 



{p.kf = — (mi + (1 - vl) V - 2ml (1 - «M • (35) 



1 

7~4 V""7T 1 "7TJ ~q TT\~ "TT/~q 

In vacuum, this expression reduces to which is also the case for a pion with zero momentum in 

the rest frame of dense matter. Note that the matrix element vanishes if either of the photons become 
soft. This helps to tame potential divergences from the Bose enhancement factors. The emissivity 
calculation proceeds similarly to that of the electroweak decay, where we had ignored the electron mass 
in the kinematics. Rescaling energies by the temperature as x = ou q /T, y = u p /T, and if) = m n /T, we 
obtain 

A\T 7 /VF* x ftp 2 fl-vl\ 2 \ 2 pnax 1 



e " ~ 2? V 2 I dx e*-l{vl [ vl) X ) J ymin dy 1 + e-*- e~v - e«- " (36) 
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For the restricted range of y and x above, the exponential factors in the y integral can be dropped, 
which is equivalent to neglecting Bose enhancement factors. Any possible divergence from this factor is 
tamed by vanishing matrix elements and phase space. The y integral is then easily performed yielding 

^ = o7 \ , / dx x e x Jx 2 -tfj 2 (^ - { y JL ) x )- ( 37 ) 



2 7 7T 3 v 3 Ji> V V V 2 . \ V 2 j / 

With the substitution x — if) cosh 9, the above integral can be expressed as a combination of modified 
Bessel functions K v (if>), in the limiting case of m n /T ^> 1. As their argument is very large, we can 
utilize the asymptotic forms of K v to obtain the result 

^ 64tt v*\2ir) ~ 64vr " W 

For purposes of numerical estimation, we note that for the typical central densities of (5-10)no, where 
no is the nuclear saturation density, the value of the strong coupling g is still much larger than e, so 
that 9c fl is small. The emissivity can then be expressed as 

= 5.1 x 10 40 m 4 (m 10 T 9 ) 3/2 /if 2 e ~ erg cm" 3 s" 1 . (39) 

Since thermal pions have an extremely small number density at keV temperatures, the photon emis- 
sivity is negligible despite the large cross-section for this process. 

It is amusing to estimate the mean free path of photons due to the inverse process 7 + 7 — > tt°. For 
the purpose of illustration only, we consider a thermal distribution of photons. The thermally averaged 
mean free path may be written as 

(A) = v^t) , (40) 



where the photon velocity is [39 



e 2 /i 2 



■/VT+Z, ^ = T^T^« 1 ■ (41) 



Setting vxj « c and utilizing the thermal average of the inverse rate 



d 3 p 



nl a;. 



(r) J (2tt) 3 t{uj p ) 
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(42) 



where 1/t(u) p ) is the typical inverse lifetime of a photon interacting with other thermal photons to 
form a 7r° excitation, we find 

<A> = ^r^ • (43) 

This result clearly displays the dependences on the vacuum decay rates and on the baryon density 
(through f n and in A^) and temperature (through and n n ). Inserting the expression for the 
number density of thermal photons, = 2£(3)T 3 /7r 2 , into Eq. (|4^), we obtain 

(A) = 1.89 x 10~ 7 £ m Ts /2 m^ ,2 e 11 ^ cm . (44) 

For the temperatures of interest, this number is exponentially large. Thus, photons in the CFL phase 
are extremely unlikely to recombine into a tt°. 

Unlike photon-photon interactions, the processes that contribute dominantly to the mean free path 
are Compton scatterings off charged mesons in the CFL phase. Since photons stem from the decay of 
thermal pions, most of the photons would have energies well below mcB, since T <C viigb < A. Thus, 
the Compton cross-section may be taken to be adequately represented by the Thomson scattering limit: 

^~°t(—) 2 ~^t I c4 = |(^) 2 = 66.5fm 2 . (45) 
\mGB/ 400 3 \m e c 2 J 

The mean free path from scattering off charged pions can be estimated as 

I ~ — - — ~ 2 x 10 40 km . (46) 

Since the density of scatterers, charged Goldstone bosons, is exponentially small, the mean free path 
is very large. Hence, at low temperatures, the CFL phase can be regarded as being transparent to 
photons [JT^] even in the presence of Goldstone bosons. 

The mean free path of a photon is essentially the size of the CFL quark core of the neutron star, 
since photons equilibriate rapidly by interacting with the surrounding ordinary matter. For normal 
nuclear matter at a density % = 3n , with n = 0.16 fm~ 3 , a similar estimate as above yields 

I' ~ -s- 10~ 8 cm . (47) 

erf n B 



We conclude that both photon emissivities and opacities in the CFL phase are negligible at the 
temperatures relevant for the long-term cooling of a neutron star. 
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IV. SPECIFIC HEAT OF CFL MATTER 



In order to determine the cooling history of a neutron star, we also need the specific heat of CFL 
matter. The specific heat of ordinary quark matter is dominated by the quark contribution 

Cq = V 7T 2 n/^ =2.44 xl0 20 T 9 f ^V /3 erg cm^K- 1 , (48) 

f=u,d,s V J KU ° J 

where for simplicity all quarks have been taken as massless and thus fif is the same for all quark flavors. 
The gluon contribution is 

4tt 2 

c = N g T 3 = 3N g x 10 13 Tg ergcnT 3 K" 1 , (49) 

15 

with N g = N 2 — 1=8. The photon contribution is identical to the gluon contribution with JV 7 = 1. 

In ordinary quark matter, there is a non-zero electron chemical potential. As a result, electrons 
provide a significant contribution to the specific heat: 

f k BT\ _ n K ~ in 2o rp fY e n B \ 2/3 3 ! 



c e = 7r 2 n e = 0.56 x 10 2U T 9 erg cm" 3 K . (50) 

V pL e J V n J 

In CFL quark matter, the contributions from quarks are suppressed by exp(— A/lksT)), where A 
is the gap in the appropriate channel. In addition, Y e 1 in the CFL phase. Note, however, that at 
T 9 = l,c e > c g + c 7 till Y e = 10- g . 

The specific heat of CFL matter receives significant contributions from the Goldstone modes. The 
contribution of massive thermal bosons (that obey the dispersion relation u 2 = v 2 k 2 + m 2 ) in the CFL 
phase is modified by the factor l/t> 3 relative to the case with iv = 1. For example, for pions in the 
non-degenerate limit, we have 

3 / 5 T \ 116m 10 

c GB {m ^ 0) = -n n 1 + -— + •••= 7.12 x 10 15 (m 10 T 9 ) 3/2 e erg cm' 3 KT 1 . (51) 

The contribution of a Goldstone mode with dispersion relation ujgb = v w k is also modified by a 
factor 1/f 3 over that of a massless mode with cu = k. With v n = l/y/3, we find 

6\/^n 2 

c GB = — T 3 = 7.8 x 10 13 Tg erg cm" 3 K" 1 . (52) 

15 
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At temperatures on the order of several MeV, there are 10 almost massless Goldstone modes. When the 
temperature drops below 1 MeV, the specific heat of CFL matter is dominated by the exactly massless 
U(1)b Goldstone boson. If CFL matter is kaon condensed, then there are two massless Goldstone 
bosons. 



V. PRINCIPAL FINDINGS 

We have studied neutrino emission from the decay and scattering processes of Goldstone modes 
in high density quark matter in the superconducting color-flavor locked (CFL) phase. Such a phase 
might occur in the core of a neutron star whose mass is near its maximum allowed value. 

The neutrino emissivities from the new processes that we have identified scale as 

K ± _> i± + Vl . €k± „ (sin 2 dc ai f 2 m 2) E 2 r± Uk± 

x i 1+2 (l ± M) + y(| ± 



A e k± ) v 2 K m 2 e \ f w E K ± 
■K°(r},r}') -> v + v : e vV ~ (G^m^Sf 2 ) T n n 




2 



+ 0^0 + z/ + z>: e^^-^T 15 , (53) 

where 5 f n = fx — fs with f s = f n /3, and and uk are the number densities of pions and kaons, 
respectively. Since the masses of these Goldstone bosons are expected to be of order 10 MeV and 
T < 0.1 MeV, these Goldstone bosons are non-degenerate and obey Boltzmann statistics. Hence, their 
number densities scale as (mT) 3 / 2 exp(— m/T). This is in contrast to quarks, which are degenerate since 
Hb/T 3> 1. However, emissivities from gapped quarks are suppressed by a factor exp(— A/T), with the 
gap A being of order 100 MeV in the CFL phase. The hierarchy of physical scales, fiB > A > itlgb > T 
(for the most part, we have taken /i e = in the CFL phase) thus assures that the neutrino emissivities 
from massive Goldstone modes dominate over those from gapped quarks. 

The unique feature of the tt° — > v + v process is that it requires Lorentz symmetry breaking, which 
manifests itself in the difference of the spatial and temporal pion decay constants, fs and fx- This 
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emissivity is thus proportional to 5f%. In vacuum, helicity conservation forbids such a decay, but in 
matter where boost invariance is lost, this process is allowed for all non-zero pion momenta. 

The emissivities from the massive Goldstone modes are expressible in terms of simple products 
interaction strength x relevant energies x the number density of thermal bosons, mainly because 
bosons are non-degenerate. Such simplicity is lost when these modes exist in partially degenerate or 
degenerate conditions because of Bose enhancement effects. It must be stressed that the emissivities 
from the massive Goldstone modes, even though they dominate those from gapped quarks, remain 
rather small, chiefly because of their very small number densities at the temperatures of relevance to 
long-term cooling. For example, at Tg = 1, n n ~ 3.5 x 10 -58 fm~ 3 . 

Neutrino emission from the scattering of massless Goldstone modes (numbering one in the CFL 
phase and two in the kaon condensed phase) is not exponentially penalized. However, this process 
appears only at 0{p A ) in a chiral effective Lagrangian, which leads to a numerically small emissivity. 
Nevertheless, it is likely to be the principal source of neutrino emission from the CFL phase. 

We have also calculated the photon emissivity from the CFL phase. The anomalous decay process 
7T° ^7 + 7 yields an intuitive expression for the emissivity, = (A^m A / 64n) n n , when m w ^> T. 
This emissivity may be readily understood as the product of the energy per unit time carried by the 
photons times the number density of thermal pions from which they are emitted. Numerically, the 
emissivity is exponentially small, not because the cross-section for the anomalous process is small, but 
because, at T 9 = 1, the thermal pion density is exceedingly small. Our estimation of the photon mean 
free path, due to Compton scattering in the Thomson limit, shows that the CFL phase is transparent 
to photons even in the presence of Goldstone bosons. However, photons are thermalized quickly in the 
high density normal baryonic matter that surrounds the CFL phase. 

Since electrons are strictly absent in bulk CFL matter, its specific heat per unit volume cy is 
dominated by the massless U(l) mode associated with superfluidity (two such modes exist for the 
kaon condensed phase). This bosonic excitation (0) is relatively easier to excite than all other massive 
Goldstone modes, which are exponentially suppressed by thermal occupation factors. Thermal gradi- 
ents arising in the CFL phase are rapidly reduced due to the large thermal conductivity, which is a 
consequence of small cross-sections for 00 scattering as well as their low number density. As expected, 
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the transport properties of the CFL phase are controlled essentially by its superfTuid nature. 



VI. IMPLICATIONS FOR THE COOLING OF COMPACT STARS 

Our results for the neutrino emissivity in the CFL phase lead naturally to the issue of their efficiency 
in cooling the interior of the star. The thermal evolution is determined by the equations of radiative 
transport and energy balance which may be cast to read as (48 



where T is the temperature, r is the radial coordinate (general relativistic corrections have been 
suppressed for simplicity), K is the thermal conductivity, and cy is the total specific heat per unit 
volume. The total emissivity due to photons and neutrinos is denoted by e and possible internal heating 
sources are contained in H. We can estimate the cooling timescales At of the CFL phase in isolation 
with the expression 

At = - f Tf dT . (55) 
JTi e 

In a strictly electron-free CFL phase, the largest contributions to the specific heat comes from those of 
gluons, photons, and the massless Goldstone modes, since those of the massive Goldstone modes are 
severely suppressed. For the same reason, the emissivity is dominated by the scattering of massless 
Goldstone modes. Using cy ~ 3.5 x 10 14 Tg erg cm -3 K -1 and e ~ 10~ n x T 9 15 erg cm -3 s _1 from Sees. 
II and III, we find that At ~ 10 26 T 9 _11 /i^ 00 y. This time scale is extremely long! 

Note that, were the massive Goldstone bosons to operate in isolation, the cooling time would be 
vastly different. Using cy ~ 7.2 x 10 15 (m 10 T 9 ) 3 / 2 exp(— 116m 10 /T 9 ) ergcm _3 K _1 and e ~ 1.2 x 
T£ s fi 2 W0 mj {m w T 9 ) 3 / 2 exp(-116m 10 /T 9 ) erg cm" 3 s" 1 , we obtain At = 6 x 1(T 4 x (T 9 (i) - T 9 (/)) 
Mioo m w s - I n this case, c v /e becomes independent of T, which leads to a cooling rate that is linear in 
T, and also exceedingly rapid. 

The CFL phase is characterized by a very large thermal conductivity due to the smallness of both 
the Goldstone boson number densities and cross sections in CFL matter. The temperature gradient in 
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the CFL core will thus be very small, and the CFL core will be nearly isothermal, with its temperature 
pinned to that of the surrounding hadronic matter (this is easily verified by seeking separable solutions 
to Eq. (0))- In the realistic case in which a neutron star has a small CFL core, the thermal content 
and the emissivity of the CFL core will be negligible compared to those of the surrounding hadronic 
matter .0 The picture that then emerges is that such a star will cool nearly identically to a star without 
a CFL core. The presence of a CFL phase in dense matter will apparently leave no observable trace 
in the long-term cooling history of the star. 

The interface between electron-rich nuclear matter and electron- deficient CFL quark matter was 



analyzed in Ref. |33| with the conclusion that a thin charged layer (~ 10 fm) separates the two bulk 
phases if the surface energy a s is large (dimensional estimates place a s ~ 400 MeV, which is too costly 
for a mixed phase). In this case, electrons can leak in to the CFL phase. Owing to the partially 
inverted mass spectrum of Goldstone bosons in the CFL phase (in particular, < rn n ), electrons 
can decay into kaons through electroweak interactions near the CFL side of the interface, resulting in 
neutrino emission and the formation of a K~ condensate. This causes the proton concentration on the 
nuclear side of the interface to rise substantially in order to maintain charge neutrality. Consequently, 



3 In baryonic matter, the simplest possible v emitting processes are the direct Urea processes f\ + i — > 
/2 + vi , f'2 — > fi + t + Vg, where /1 and /2 are baryons and £ is either an electron or a muon. These processes 
can occur whenever momentum conservation is satisfied among fi, /2 and £. If the unsuppressed direct Urea 
process for any component occurs, a neutron star will rapidly cool because of large energy losses due to 
neutrino emission: the star's interior temperature T will drop below 10 9 K in minutes and reach 10 K in 
about a hundred years. This is the so-called rapid cooling paradigm |4S,5C]. If no direct Urea processes are 



allowed, or they are all suppressed due to baryon superfluidity (with gaps of order 1 MeV), cooling instead 
proceeds through the significantly less rapid modified Urea process in which an additional fermion enables 
momentum conservation. This situation could occur if no hyperons are present, or the nuclear symmetry 
energy has a weak density dependence [49, |0| . Comparatively less rapid, but faster than that due to modified 



Urea cooling occurs in the presence of Bose condensates 
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there is the possibility that, within this very thin layer, the direct Urea process n — > p + e _ + v would 
occur as long as the proton fraction is of order 11 — 14% for which momentum conservation between 
the participating fermions becomes possible |[49|| . Inasmuch as this thin interface is embedded within 



the star, the occurence of a such a rapid cooling process would be masked by cooling from the much 
thicker exterior layers in which similar processes are also possible. From an observational standpoint, 
such cooling would resemble cooling from normal stars. 



In Ref. |33|, the phase structure of baryon and CFL quark matter was also studied as a function 
of the magnitude of a s . For a sufficiently small a s , a homogeneous mixed phase with domains of CFL 
quark matter and nuclear matter was found to be preferred over a wide range of density. In this case, 
neutrino emission would be dominated by the nuclear bubbles which can cool via the rapid processes 
characteristic of the nuclear phase (see footnote 3). 



Recently, the cooling of a self-bound quark star with a bare surface has been studied in Ref. |5T 
by considering both the CFL and 2SC phases. If the entire star is in the CFL phase, which permits no 
electrons, cooling occurs chiefly through processes that involve the Goldstone modes. As noted above, 
the cooling times for this case are too long to be observable. If the CFL phase is surrounded by the 
2SC phase at lower densities, cooling is dominated by neutrino emission from direct Urea processes 
involving unpaired quarks in the 2SC phase (which allows electrons to exist) and by the annihilation 
of e + e~ pairs produced by an intense electrical field which binds electrons to the surface of quark 
matter. This is essentially the Schwinger mechanism in the presence of a Fermi sea of electrons, and 
can occur only at finite temperature [52]. Beginning at a temperature of ~ 10 11 K (~ 10 MeV) down to 



~ 10 8 K, the temperature versus age curve for the case with the 2SC phase in which gaps are of order 
100 MeV differs little from that of the normal phase, chiefly because pairing reduces neutrino emiision 
and the specific heat by similar amounts. Substantial differences from the cooling of the normal phase 
occur, however, if quarks not participating in 2SC pairing also pair through residual interactions, but 
with gaps of order 1 MeV. In this case, the hard X-ray spectrum with a mean energy of ~ 10 2 keV is 
proffered as an observational signature of bare quark stars. 

In this work, we have focused on the long-term cooling epoch when most neutrinos have left the 
star and temperatures are in the range of hundreds of keV to tens of eV. In a proto-neutron star, 
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neutrinos are trapped and temperatures range from 20-40 MeV. The CFL phase is likely to appear 
during the PNS evolution. During this stage, mcB/T < 1 and hence the Goldstone modes are in 
the partially degenerate or degenerate regime. As mcB/T < A/T, both neutral and charged current 
neutrino interactions with the Goldstone bosons offer new sources of opacity. Some examples of possible 
processes are z/(z/) + 7r ± — > z/(z/) + 7r ± , z/+7r~ — > e~ + tt°, v+tt + — > e + + ir° and u(u) + 7T° — > e T + ir ± . Note 
that with fast neutrinos, the processes v — > 7r° + z/ and v — > 7r + + e~ also become allowed and contribute 
to the opacity. The resultant neutrino opacities are expected to be significantly larger than those from 
scattering and absorption on gapped quarks, since temperatures are high enough to admit large pion 
densities. For example, in the range 20 < T/MeV < 40, the number densities of a given species of 
a Goldstone mode whose mass lies in the range < m GB /MeV < 10 are ugb ~ 10~ 3 — 10~ 4 fm~ 3 , 
substantially more than those encountered in the long-term cooling epoch. 



NOTES: (a) Concurrently and independently of our work, Reddy, Sadzikowsky, and Tachibana |53 
have performed a detailed study of //-opacities and emissivities in the CFL phase at tens of MeV tem- 
peratures encountered in the evolution of a proto-neutron star. Their quantitative results underscore 
the important role of Goldstone boson excitations in the CFL phase and provide the groundwork for 
calculations of observable neutrino luminosities to be performed, (b) Shovkovy and Ellis have 



performed a calculation of the thermal conductivity in the CFL phase at temperatures of relevance to 
long-term cooling. They find that the dominant contribution to the conductivity comes from photons 
and Goldstone bosons. This result is a necessary ingredient in quantitative calculations of a compact 



star's long-term cooling, (c) Page and Usov [51] have calculated the thermal evolution and light curves 
of young bare (self-bound) quark stars considering both the 2SC and CFL phases. They find that the 
energy gap of superconducting quark matter may be estimated from the light curves if it is in the range 
~ 0.5 MeV to a few MeV. 

We would like to thank P. Bedaque and J.M. Lattimer for useful discussions. We are grateful to 
K. Rajagopal for his comment concerning the thermal conductivity of the CFL core and to S. Reddy 
for helpful remarks with regard to emissivities involving kaons. This work was supported in part by 
the US DOE grant DE-FG-88ER40388. 
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